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ABSTRACT 

An inclined solar still, combining the process of 
desalination and hot-water production has been developed. 
Theoretical analysis has been given. Tests have been carried 
out for a summer month and performance evaluated. The yield 
of the inclined still has been compared with that of a standard 
basin-type unit under similar condition. 

The essential feature of the solar still developed in 
the present work consists of a thin film of water flooring 
continuously over the glass cover, which keeps it cool resulting 
in high yield of distillate. Heat gained by the water as a 
result of the heat transfer from the saline water in the basin 
to the glass cov^r raises its temperature as it flows over the 
cover. Thus, hot-water is obtained as a by-product in this 
process. The entire unit is kept at optimum tilt from the 
horizontal to receive maximum amount of incoming solar radiation 
on the basin surface. 

A theoretical analysis of the inclined solar still has 
been carried out and the results compared with those obtained 
experimentally. 

It is concluded that the inclined solar still with 
cooling-water system is an efficient device for the 
purpose of desalination. Such units may be installed over 
large areas near sea-coasts for large-scale production of 
freshwater. 
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IITTRODUCTION 


1.1 WATER SCARCITY, A (E.0BAL PROBLj^'l 

Lack of adequate quantities of fresh water has 
become a serious problem in several countries around, the 
world. Large areas of the earth are arid, due to over- 
whelming sunshine and lack of rains. In most of these 
regions, ground water is available, but with a high 
percentage of mineral content, which makes it unfit for 
drinking and agricultural purposes. Small villages, 
situated near sea-coasts, too, face shortage of fresh water 
during peak periods of summer w-hen fresh-\ 7 ater rivers and 
lakes nearby become dry, rendering serious drinking water 
problem to the population. ’Demineralization process’ to 
produce desalinated water from saline or brackish water is 
not possible to carry out in most of the remote areas, due 
to non-availability of electricity, necessary for pumping 
processes and also due to high cost of installation of the 
demineralization plant. 

Thus, in most of these arid zones or regions near 
a sea, where sunshine is abundant, coriversion of saline or 
brackish water (also available in large quantities) to 
fresh water by means of solar energy seems to be the only 


possible solution today 



2 


1.2 PRBJCIPLS OP SOLAR DESIJjJE AT1(M 

The principle of conversion of saline or brackish 
water to fresh water by means of solar energy can be best 
■understood by referring to the simplest device, designated 
as the “basin-type solar still."' This device is a repro- 
duction, oa a very small scale of the hydrologic cycle that 
occurs in nature, namely, heating of lakes and oceans by the 
sun's energy causing evaporation, the transport of water- 
vapour to cooler regions by convective winds and the subse- 
quent condensation of the vapour causing rain and snow. As 
shown schematically in Pig.l,l,it consists of a blackened 
basin containing saline water at a shallow depth, ever which 
is an air-tight glass cover that completely encloses the 
space above the basin. Solar energy, passing through the 
transparent cover, is absorbed by the 'water and the basin, 
warming the water causing evaporation to take place, A 
temperature difference exists between the hot saline-water 
surface and the relatively cool glass cover, resulting in 
convection currents to set up in the 8.ir, trapped inside 
the enclosure. These currents bring the humid air in 
contact with the glass cover, causing condensation of the 
vapour on the underside of the covei’. The condensate 
slides down the slope, flows along the distillate channel c and 
is collected in a container. 



GLASS COVER 
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1.3 ADVMTAGfiS SOLAR DESALIRATIOR 

Some of the important factors fa.vouring solar 

desalination are as following s 

(1) In arid regions where there is no other source to 
obtain drinking water, solar desalination is the 
only solution, 

(2) Solar energy is used as a substitute for fuel in 
solar desalination process to heat the saline water 
to form vapours which later condense to produce 
fresh x^ater. In view of the present scarcity and 
hi^ cost of fuels like coal, oil and gas all over 
the world, use of solar energy, which is free, 
inexhaustible and non-polluting in nature, is an 
attractive idea for desalination purpose. 

(3) Solar desalination is most effective during the 
peak periods of sunshine, when the demand for fresh- 
water is also every high. Also, during rainy season, 
when solar stills are not very effective, they can 
be used to collect rain water as fresh water. Thus, 
the donand and supply are effectively geared in this 
desalination process, 

(4) Unlike the problem of installation of flat-plate 
collectors in populated areas where space problen is 
acute, solar stills in most cases, are installed 
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over vast areas in the arid region, which are other- 
wise lying barren. Production of fresh water in 
these arid zones may facilitate in utilizing large 
areas of land In future. 

1.4 REVIEW OF PREVIOUS WORK 

Historically, instances of application of solar 
energy 'to fresh water production have been recorded 
throughout the world and at various times during the past 
century. The unit which has been built and used exten- 
sively in the past and is still the most popular device 
for desalination purpose today, is the simple basin-type 
solar still, already discussed. The largest basin-type 

solar still was built in the year 1872 near Las Salinas in 

2 

northern Chile. This unit had a basin area of 4700 m and 
it produced 22700 litres of pure water per day to the needs 
of a nitrate mining ccmmunity for many years. 

During the World War II and the decade immediately 
follovjing, sustained drou^t conditions in many parts of 
the world brought water suuply into sharp focus. This 
i-esulted in the development of inflated plastic stills by 
Maria Telkes^. The unit consisted of an inflatable plastic- 
envelope containing a felt-pad which was to be saturated 
with sea-water. Solar energy, striking the felt-pad, would 
produce vapour which would condense on the inside of the 
plastic-envelope and drip into the bottle at the bottom of 
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the assanbly. Althou^^, plastic stills are relatively 
cheap j their yield is low due to the formation of drop- 
X'j’ise condensation under the plastic roof, causing hi^ 
reflection losses of incoming solar radiation. The effi- 
ciency of such stills, is less than 40 percent. 

Tilted solar stills using porous pads as evapora- 
ting surfaces have also "been used widely for many years. 

The upper end of the pad is immersed in a trou^ containing 
salt water, which spreads over the cotton by capillary 
action. Evaporation and condensation take place as usual. 
Although, the unit is cheap and easy to fabricate, yield is 
quite low due to low temperature of the evaporating surface. 
A maximum efficiency of 50 per cent has been reported for 
such stills. 

A solar still with a low thermal inertia has been 
2 

developed by Ssulmayer . It is essentially a basin-type 
still with a floating, black absorber sheet to heat a thin 
layer of water (1.6 - 6.3 mm thick) on the top of the sheet, 
causing rapid evaporation. Since the base of the still is 
kept cooler, heat loss throu^ conduction to the ground is 
reduced. The rate of evaporation (using a woven black shade 
cloth as the floating absorber^ has be^iui reported to be 
D. 3 O 8 litres /m^-hr. . 

Efforts have been made to increase the temperature 

of the basin-water by means of concentrating devices. A 

■ 3 

'mini solar-still’ designed by Muthuveerappan, et al. , 
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incorporates mirror-type reflectors to concentrate solar 
energy on the tray. A 50-60 per cent increase in the yield 
of distilled water and a mean efficiency of 56 per cent 
have 'been reported, Hox^ever, the extra cost of reflectors 
seems to outweigh the effect of increase in the yield. 

4 

Loho, et al. have tested a multi-effect, basin- 
type solar still, which uses the upward heat fluxes due to 
convection, radiation and evaporation to heat a secondary 
basin consisting of an inclined- stepped glass tray. The 
upward heat-fluxes from the secondary basin are used to 
heat water contained in a third basin, similar in construc- 
tion to the previous one. The condensate frcm each effect 
is collected in separate bottles, Althou^, an increase in 
yield of 40 - 50 per cent over the standard basin- type unit 
has been reported, the unit suffers from constructional 
difficulties, leakages of both saline-water and -water-vapour 
throu^ the linings and lack of long-term durability. 

A second type of regenerative unit has been developed 

5 

by Akhtamov, et al. . This unit incorporates the reuse of 
upward heat fluxes to heat a layer of saline-water which is 
enclosed between two glass-covers and flows upwards due to 
hydrostatic pressure exerted by the saline-water contained 
in an overhead tank. The saline-water absorbs heat and is 
allowed to be fed ijito the basin to ecmponsate for the 
evaporation of saline water. This process results in the 
preheating of the saline feed-water to' the basin. An 
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efficiency of 60 per cent has been reported by the authors. 
However, since the rate of discharge of hot saline water to 
the basin is very small, the problon of cooling the cover 
assembly to obtain higher yield, yet remains to be solved. 

In view of the large area required to produce a 
given quantity of water by means of the basin-type stills, 
investigators have tried vertical type, multi-effect solar 
stills. One such unit developed by Dunkle^ consists of 
vertical parallel plates within an enclosure in which hot 
saline water is made to flow down the sm-face of a plate 
causing evaporation. Condensation take splace along the 
surface of an adjacent vertical plate, vrhich is kept cool 
by allowing cold saline water to run dom along the other 
surface of it. A series of such evaporators and condensers 
are installed to obtain higher yield. However, this increase 
in distiller performance is yet to be justified in terms of 
the large capital cost on the total plant economy. 

1,5 THE PRISM T WORK 

In most of the work on solar desalination, as 
reviewed above, the conversion efficiency has not been 
more than 60 per cent. The present effort is directed 
towards developing a more efficient still which acts as 
an auxiliary hot-water device as well. The main feature of 
the solar still developed in the present work consists of 
a box- type unit in vrhich a thin film of water flows 
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continuously over the glass cover and keeps it cccl, 
resulting in high yield of distillate. The saliiie water 
is stored in an overhead tank and fed continuously to the 
still through a perforated tube. Upward heat fluxes due 
to convection, radiation and evaporation get transferred to 
the thin film of water, whose temperature increases as it 
flows along the glass. Thus, in addition to high output 
of desalinated water due to increased evaporation rate, 
hot water, also, is obtained continuously throughout the day, 
which may be either stored or used to heat some auxiliary 
devices. The entire unit is kept at optimum tilt frcxn the 
horizontal to receive the maximum amount of incoming solar 
radiation on the basin surface. 

Chapter-2 describes the design and fabrication 
details of both the inclined and the basin-type stills for a 
comparative study. It also describes the instrumentation 
involved in the experiments and the modes of tests conducted. 

In Chapter-3, heat and mass transfer analysis of 
the inclined still has been carried out. 

The resifLts obtained both from theoretical anal^^sis 
and experiments carried out are discussed in Chapter-4. The 
hourly and cumulative yields of the inclined still are 
compared with that of the standard single roof, basin-type 
unit for a summer month. 



CHAPT-3R ~ 2 


THE EXPERmaTTiOL SET-UP 

Tile experimental set-up that was used for testing 
the performance of the inclined solar still and for compa- 
ring it with the oasin-type unit is shown in Pig. 2.1. 

2.1.1 DESCRIPTION OP THE niCLBiED SOLlR STILL 

The constructional details of various important 
parts of the inclined solar still are shown in Pig. 2.2 and 
discussed below. 

The saline-water basin consists of 60 cms long 
aluminium angles (1cm x 1cm, 0.5 mm thick) placed in 
adjacent rows and fixed to the bottom of a rectangular 
box (60cms x 60cns x 10,6 cms) made of 22 gauge (0.85344 mm 
mean thickness) G.I. sheet. The channelled basin, thus, can 
contain about 3.6 litres of saline water when filled to a 
depth of 1 cm. To prevent leakage of water frcm one channel 
to the other throu^ the sides, a rubber gasket, 1 mm thick, 
surrounds the periphery of the basin. The basin is painted 
with black paint mixed with black mesh powder, and later 
made dull by rubbing with emery paper to make it a good 
solar energy absorber. 




FIG.21 SCHEMATIC DIAGRAM OF THE EXPERlMENTAi. SF T UF 
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Common window ^ass (60cms s 60cms, 5.5 mm thick) 
is used for the top covers of the still. Glass transmits 
incoming solar radiation in the range of 0.3 to 3»0 microns, 
hut is opaque to radiation emitted by the hot absorber plate 
and the water surface, beyond 4.8 microns. Hence, the still 
serves as a hot-box. The glass cover rests • on aluminixM 
angles, fixed along the periphery of the box. It is sealed 
to the box by means of rubber gaskets and adhesive along 
its four sides to make it air-tight, so that no water-vapours 
may leak outside throu^i the edges. A gap of 5 cmsis maintained 
between the basin water surface and the glass cover. This 
distance of 5 cms was chosen to reduce the shading of the 
basin water surface by the lateral walls of the solar still 
and also to decrease the fabrication cost. The top glass 
cover, similarly, rests ligh'fly on aluminiiam angles, provided 
with rubber gaskets, along the periphery of the box, but it 
is removeable. The top glass cover is periodically treated 
xfith dilute hydrochloric acid to rancve the formation of 
dust coatings. 

To prevent heat losses from the bottom and sides of 
the still, the unit is enclosed in a rectangular box 
(64cms X 64cms x 16.6cms) made of 22 gauge 6.1. sheet, and 
the space in betvreen is filled with 6 cms thick ^ass wool 
at the bottom and 2 cms thick on the sices. 
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2.1.2 IHj5 BASUT-TTPE SOLAE SXILl 


The c on struct! on al details of the hasin-type still are 
shown in Fig. 2.3. 

It consists of a box with base dimensions of SOcms z 
60 cms, front and rear wall heights of 6 cuts and 18.5 cmSj 
respectively. The still has an inclined roof on which rests 
a 5-5 mm thick glass cover, sloping at an angle of 12*^ with 
the horizontal to prevent the condensate from dropping into 
the basin. The glass cover rests on rubber linings and is 
secured properly in its place by means of aluminium angles, 
fixed with screws along the periphery of the still. The 
condensate slides dom the slope of the glass cover, flows in 
a U-shaped channel from where it is collected in a plastic 
bottle. Saline water is filled upto a small depth (about 1 cm) 
in the basin for maintaining a low thermal inertia. The basin 
surface is painted with black paint mixed with black mesh powder 
for better absorption of solar radiation. 

The basin- type still is also enclosed in a box with base 
dimensions of 64cms x 64cms, front and rear wall heights of 
12 cms and 24.5 cms, respectively. Both the inner and the 
outer boxes are made of 22 gauge G-.I. sheet and the space in 
between is filled with 6 cms thick glass wool at the bottom 
and 2 cms thick on the sides to prevent heat losses to the 
surroundings. 
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2.1.3 PIPINGS, FITTINGS AND STORAGE TMK 

The arrangement of the pipings, valves and other 
connections for both the stills is shorn in Fig. 2.1. Polythene 
pipes with plastic valves and connections have been used for 
ease of installation and low cost, as well. For the inclined 
still, the condensate is collected in a U-shaped channel, 
running along the underside of the glass cover at its edge 
as shown in Fig. 2. 2. Also, a 1 cm dia. brass tube with uniformly 
spaced holes of 2 mm dia. drilled along its length is used to 
feed the cold saline water to flow over the glass plate. The 
hot-water obtained as a by-product, flows out through the exit 
pipe into a trough placed on the ground, from where it is 
drained away. 

A storage tank of capacity 200 litres, is provided 
to feed the saline water to both the stills. It has two 
separate valves connected to the outlet pipes of the tank, 
providing saline wa^tertoeach basin and cold saline water to 
cool the bottom ^ass cover of the inclined still. A graduated 
water-level indicator made of glass is fitted to the storage 
tank. The tank is insulated with 6 cms thick glass-wool along 
its sides to prevent the stored saline-water from being heated 
up. 
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2.2 lUSTOTiENTATION 

2.2.1 TOm SOLjffi RADIATION MEABUREMENT 

The total radiation incident on the basin of the 
inclined still consists of both the direct and diffuse radia- 
tion. An Eppley Pyranometer has been used to measure the 
total radiation falling on the inclined surface. The instru- 
ment was kept facing south and tilted at the angle of incli- 
nation of the still, and the output recorded on an hourly 
interval by the Eppley Integrator. 

2.2.2 TEMPERATURE MEASUREiyrENT 

All t super atures, ■ except the ambient have been measured 
by 24 gauge, calibrated copper-const ant an, insulated thermo- 
couples, and the output recorded directly on the graph paper 
by a 24“3unction Honey-well Temperature Recorder, having a 
least count of l^C . Fig. 2.4 shows the positions of the thermo- 
couples on the inclined still. The top ^ass cover has four 
thermocouples fixed to its upper surface by adhesive as shown. 
Similarly, four thermocouples have been fix;ed to the underside 
of the bottom glass cover to measure the average temperature of 
the condensate film. Four thermocouple wires have been soldered 
to the channelled basin and four to the bottom surface of the 
still. Two thermocouples, each immersed in the flowing water 
film were used to measure the inlet and outlet temperatures 
of the wa'fcer film. Ambient temperature was measured by a mercury 
in-^ass thermometer with a least count of 1°C, placed in shade. 




2.2.3 PLOW MEASURMMT 


A graduated, water-level indicator made of ^ass, 
fitted to the storage tank, as shown in Pig. 2.1, is used to 
note the decrease in the level of saline water over an hourly 
interval of time. Enomng the area of cross- seotion of the 
tank, the mass flow z’ate of saline water oan be found out. 

Ihe condensate collected in plastic bottles was 
measured on an hourly basis by means of a graduated cylinder 
of 500 ml capacity, having a least count of 5 ml. 

2.2.4 WIND VELOCITY MEASURS'iSNT 

Due to non-availability of a windsoope, the velocity 
of wind was measured by means of a rotating, vane- type anemo- 
meter, which was adjusted every time to take into account the 
change in the direction of the wind velocity. 
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2.5 MODES OE TEST 

The performance of the inclined still was tested 
■under t-wo modes of operations. 

In the first mode, the top glass cover (-which is 
removeahle) was placed over the still to prevent cooling of 
the film of flowing water due to its surface leing exposed 
to the ambient, thus retaining the hot water which often 
reached more than 48 °G during the peak hours of sunshine. 

This arrangonent , however, decreased the amount of incoming 
radiation due to the addition of the top ^ass cover, resulting 
in lower yield. 

The second mode of test -was conducted with the top 
glass cover removed. In this case, yield of distillate was 
more, but the temperature of the hot water was substantially 
reduced due to the rapid heat loss from the water film to 


the ambient. 



CHAPgjjR - 3 


HEAT AND MASS TRAHSEER ANALYSIS OP THE 
IHCLIHED SOLAR STILL 

3.1 HEAT BALANCE CF THE STILL 

Assijming steady-state operation of tlie still, th.e 
following equations can be written : 

Heat balance for basin water : 

^a = Rc ^ ^r 


Heat balance for the bottom ^ass cover : 


lo + “Ir + le 


m 

A 


g 


c (T 
p w 


- T ) 
o ^i 


+ q 




(3.2) 


where, 

^a 

^c 

^e 


incident solar energy fltnc absorbed by the 
basin-water, 

heat flux due to natural convection from the 
basin-water surface to the ^ass cover, 
heat flux due to radiation frcaa the basin-water 
surface to the glass cover, 

heat flux due to evaporation atid condensaticn 
from the basin-water to the ^ass cover, 



22 


heat flux due to hack losses irora the basin 
to the surroundings, 

heat flux due to upward losses from the ^ass 
cover to the surroundings, 

mass flow rate of saline water over the glass 
c over , 

area of the glass cover, 

- inlet and outlet temperatures of the saline 
water flowing over the glass cover, 
respectively. 

The various heat fluxes are evaluated in the following 
sections : 

3.2 HEAT FLUX DUE TO NATURAL CONVECTION 

The process of heat transfer inside the still is due 
to natural convection arising as a result of the density 
difference of the fluid enclosed between the two parallel 
surfaces, the basin surface at temperature T^ and the glass 
surface at tanperature T^y where T^ > T^. Under such circums- 

tances, the heat transfer coefficient is obtained from the 

7 

anpirical relation , 

Nu = [ 2. 44+ { 2 . 23-0. 01236+0. 00034©^} - ij-] 

valid f or 10 ^ G-r Er <C^4xl0^ ' 

0® < 6 < 60° (3.3) 


% 




m 


T ,T 
^i ^0 
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wh er e , 

Nu = Uusselt num'ber 

V; 

k 

Gr = Grastiof number for heat transfer 


g 1^ 2 (T^ - T ) 



Pr = Prandtl number 



Ra^ = 11300[1 + 0.204 sin 4.50(8 - 37.8°) ] 

6 = Angle of inclination of the still with the 

horizontal 

1 = Distance between the two surfaces at tempera- 
tures and T2» respectively 

g = Acceleration due to gravity 

= Coefficient of volume- exp ansi on of the air- 
water vapour mixture in the still 

h = Convective heat- transfer coefficient 

c 

yO , p., c , and k denote the density, viscosity, specific 
heat and thermal conductivity of the mixture in the still. 

The term p + (Tn - To)» known as the ’buoyancy factor’ for heat 
transfer was originally evaluated for conditions of constant 
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composition of the conv acting fluid, and was given by 


Pt 





T, 


+ TJ/ 




However, in the still, the humidity of the air and, therefore, 
its composition changes along the flow path. The vapour 
formed in the still is lighter than the dry air component and 
therefore, increases the buoyant force and the rate of circu- 
lation of the air and associated -vapour. The modified expre- 
ssion for the buoyancy factor as derived in Appendix-1 is 
given by 


Pt'h - 



P 


M. 


( „ p 


T. 




B A 


( 3 . 4 ) 


wh er e , 


P 


An 


Pa 


2 


P 




partial pressure of saturated water-vapour 
at temperature T^. 

partial pressure of saturated water-vapour 
at temperature Tp. 

total pressure inside the still 
(assumed at 1 atmosphere). 

Molecular weight of wate3>- vapour. 

Molecular wei^t of dry air. 
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The expression (3.4) is used to determine the Grashof number 
occuring in expression (3.3). 

The heat flux due to convection is given by 

'lo = *^0 (h - ''a' 

where, the heat traxisfer coefficient is evaluated from 
equation (3*3). 

3.3 HEAT ELUZ DUE TO RADIATION 

The heat transfer due to radiation can be approximated 
as taking place between two parallel surfaces of equal areas 
enclosing a non-absorbing, non-emitting medium. The radiation 
involved Is diffuse radiation in the long wavelength range, 
for which the glass and water can be treated as grey surfaces 
(e^= e). Under such circumstances, the heat transfer due to 
radiation is expressed as 

= S' » {T^ - Tg) 

where, 

a = Stef an -Boltzmann constant 

_ Watts 

= 5.67 X 10 ^ ' ^ " "4 

m - K 

F = Shape factor 


1 
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~ emissivi'ty of the saline-water in the basin. 

Ogi = emissivity of the water film-glass glazing, 

fating the emissivities of saline-water and glazing as 0.96 
and 0,95> respectively, we obtain F = 0.9. Hence, 

= 0.9 0 ( ) (3.5) 

3.4 HEAT FLUX DUE TO EVAPORATION MD CaTDHTSATION 

In addition to convection and radiation, the bulk of 
the heat transferred upwards is due to evaporation of water 
from the basin surface and its subsequent condensation on the 
gLass cover. The evaluation of the heat transfer due to 
evaporation and condensation is carried out in three steps, 
namely, evaluation of 

(i) mass- transfer coefficient, 

(ii) rate of mass-transf er , 

and (iii) heat flux due to evaporation and condensation. 

(i) Evaluation of the mass- transfer coefficient : 

The process of heat transfer due to natural convection 
inside the still is accompanied by mass-transf er due to evapo- 
ration of water from the basin. For the case of simultaneous 
heat and mass transfer inside a fixed geometri?- and in the 
absence of any chemical reaction, it may be shown that the 
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governing equations and. the boundary conditions for heat and 
mass transfer processes are similar and hence, they can be 
treated as analogous in nature ® . The solution of heat 
transfer problem in case of natural convection can be written 
in terms of dimensionless quantities as followss 

Nu = f (Gr, Pr) 


Since the heat and mass transfer processes are analogous, 
the solution of mass-transfer problem can be written as 


Sh = f Sc) 

where, the functional relationships are the sarne^* 

The dimensionless quantities in the above equations are: 


Sh = Sherwood number 

Sc ' = Sclimidt number 


and 


Gr 


m 


A3 

= Grashof number for mass-transfer 
g (c^ - c^) 


where, h^, and D denote the mass-transfer coefficient and 

’ b _ab 

diffusion coefficient, respectively, of water-vapour in air. 
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and C2 denote the concentrations of air water-vapour mixture 
at the lower and upper surfaces of the still, respectively. 

Pjjj is known as the coefficient of densification which indicates 
the variation of density with composition for mass-transfer. 

The term p^Cc^-c^) is evaluated in a similar manner as the teimi 
p ^(1^-12) [See Appendix- 1 ] and is given hy 




T-, 


- P, 


+ [ 




] 


(5.6) 


The expression ( 3 . 6 ) is used to determine the value of Gr^. 

Hence, in equation (3.3)9 after replacing Hu by Sh, Gr by Gr^ 
and Pr by Sc and solving for the mass transfer coefficient, we get, 


[ 2. 44+ 2, 23-0.0123©+0.00034e^j*|^( 


Gr 


m 


Sc cos 0 ^ 1/3 


Ra 


4] 

(3.7) 


(ii) Evaluation of the rate of mass-transfer 

Let the two components of the gaseous mixture inside 

the still, consisting of water vapour and air be denoted by 

A and B, respectively. The molar flux of component A in the 

^ 10,11 

z - direction is given by ^ 


H 


Az “ ^ ^ 


% " "^2 


) 


I’Blm 


(3.8) 
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where P^^jn “ log-mean pressure of component B 

T = mean temperature of the mixture 
= Universal G-as Constant 


Bor normal air, water-vapour mixture at atmospheric pressure, 


the partial pressure of one component does not vary signifi- 


cantly in ccmparison to the total pressure 


12 


H enc e , 


P 




B3mi 


From Dalton’s Law of partial pressure, 


Pa + Pb = P 


Therefore, 


^Blm 


(p ” P 


A 


) + (P - Pa ^ 


1 



(3.9) 


Substituti'ng the value of p ^-j ^ in Eqn. (3*8) and denoting 
^Az obtain, 


(P 




--V (2p 




“ P, 


) /2 


(3.10) 
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Thus, h denotes the mass-transf er rate of -/rater-vapour due to 

evaporation. The value of h^ is obtained from Bqn. (3.7). 

13 14 

The properties of the mixture are evaluated ’ at the mean 
temperature, (T^ -i- '2^/2. 


(iii) Heat Flux due to evaporation and condensation : 


The heat flux due to evaporation and condensation 
is given by 




where, h^^ denotes the latent heat of condensation 
evaluate d^at the glass temperature ^ 2 ' 


A computer programme to evaluate the mass-transf er rate 
alongwith the upvrard heat fluxes due to convection, radiation 
and evaporation is included in Appendix- 2. 


3.5 HEAT mix DUE TO BACK LOSSES FRffi THE BASIN TO THE 
SURROUNDINGS 

The heat transfer due to back losses is given by, 

% = I ( h - ’ 

thermal conductivity of the glass wool insulation, 

thickness of '^ass-wool insulation at the bottom, 
temperature of the back surface of the still. 


where. 


k = 


3^b = 
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3.6 HEAT FLUX DUE TO UPWAED LOSSES EROX THE GLASS COYER 
TO THE SURROUl'IDINGS 

The rate of heat dissipation depends both on radiation 
to the sky and on convection hy air circulation. Radiation 
to the sky depends on the equivalent sky temperature, ^ 3 ]^ • 

The convective portion, is taicen as a function of wind velocity. 
The total heat flux is given by 

where, 

T = temperature of the top glass cover, 

o 

T = taaperature of the ambient, 

Si 

e = total emissivity of glass (taken as 0.95), 




The sky can be considered as blackbody at some equivalent 
sky temperature, which accounts for the fact that the 
atmosphere is not at a uniform temperature. Among the 
several expressions proposed for evaluating the sky 
tonperature, the one which is most commonly used is 
given by 


^sky 

where, and T^ 

temperature in , 


= ®a - 

denote the sky tanperature and ambient 
respectively. 



a 


Stef aXL-Jioltzmaxua constant 
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= 5.67x10 


-8 


Watts 

2 

m - ji 




outside convective heat -transi er coefficient, 
and is ^^iven by 


h; ¥ 

~T~ 


0 .0 28 Re® * ® 


where, ¥ = width of the still, and 
Re = Reynold's number 

P- 

V = Average wind velocity^ 

a* V 0 

IJ. and k denote the density, viscosity and thermal conductivity 

13 f 

of the air , evaluated at the film temperature,^— u- 


The heat flux, obtained experimentally, as shown above, 

is compared with the value calculated from Eqn. (3*2) and the 
results shown in column 8 of Table 4.1(b). 


3.6 EFFECT 17 ENBSS OF THE GLASS-COVER ASSS>iBIY AS A HEAT 

EXCHANGER 

The glass-cover assembly serves as a heat exchanger 
in which heat transfer takes place between the water-film 
enclosed by the cover ass^bly and the upward heat fluxes 
due to convection, radiation and evaporation. Neglecting the 
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small amount of heat transfer in the inlet pipes, the saline 
water enters at the outlet tonperature of the water in the 
tank, and leaves a-t a considerably higher temperature due to 
the heat gained, as it flows over the glass surface. The 
effectiveness, B of the cover-assembly as a heat exchanger can 
be defined as the ratio of the actuaa quantity of heat raaovedby the 
flowing water to the maximum possible heat transfered to the 
water. i, e. 


E 


a- c 

A p 


( T 


- T ) 
Wo w. 


Ic ^r 


3.7 ST ILL EFF IC laiCY 


A measure of the performance of a still is the conven 

sion efficiency, defined as the ratio of heat transferred in 

condensation process over any time period to the incident 

16 

solar energy over the same time period , 


3 - • 0 • 


where, 


ri 


/^e 

fh 'I" 


Ig = total incident solar radiation (direct + diffuse) 

falling on the inclined surface under consideration. 


T 


time period. 



CHAPTER - 4 


RESULTS AHL DISCUSSIONS 

4.1 EXPERIiyiEIttAL PROCEDURE 

Experiments on both the inclined and the basin-type 
stills were carried out over the month of May, 1981, and the 
data obtained for clear days of uninterrupted sunshine, only, 
are presented here. 

The storage tank was filled every morning with saline 
water at a concentration of 2.5 gms of salt dissolved per 
litre of solution. The basins of both the stills were filled 
next with saline water to a depth of about 1 cm. The saline 
water inlet valve is then adjusted to maintain a thin film of 
water flowing over the glass cover, lasting till the end of 
the day, when the tank becomes nearly empty. The glass cover 
surface was periodically treated with dilute hydrochloric acid 
to maintain the uniform flow of water and to prevent the 
formation of dry patches. Readings were taken at hourly 
intervals of time, starting from 9.00 a.m. till 5.00 p.m. 

The inclined still was tested under two modes of 
operations ; 

(i) The top glass cover was retained. This resulted in 

low yield of distillate but the temperature of hot water 
available as a by-product was high (about 44°C), 
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(ii) The top glass cover was roaoved. The yield of distillate 
was hi^ hut the hot water temperature was low (about 
39 °C). Tables 4.1? 4.2 and 4.3 present the experimental 
data and the results obtained for the two modes of tests 
carried out on typical summer days. 

4.2 ANALYSIS OF THE EXPERIMENTAL DATA iu® RESULTS 

A graphical representation of the experimental data for 
the first mode of test is shown in Eig. 4.1. The hourly total 
insolation on the inclined still is plotted alongwith the 
temperatures obtained. From the graph, it is seen that the 
maximum temperatures of the basin surface (Ij_) ? glass cover (T2) 
and the back surface (T-^^) are reached sanetime after 1 p.m., 
an hour after the period of maximum insolation. After attaining 
the maximum temperatures, the curves drop due to the decrease 
in incoming radiation and increase in heat losses from the 
basin to the ambient. 

The saline water exit temperature, shows a 

continuous rise till about 3 p.m. after which it starts 
decreasing. This may be explained as follows i 
The rate of discharge of saline water from the storage tank is 
non-linear in nature as is evident from column-3 of Table 4.1(b), 
In the morning hours, when the tank is nearly full, the mass 
flow rate is high and the upward heat fluxes are low due to low 
intensity of radiation during this period of the day. Hence, 
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Table 4.1(a) s Experimental data obtained for the inclined still 
for the first mode of test (top glass cover 
retained) ' May 26, 1981 

e = 12° 







Jable 4.1(b) t Resaits of tha test conducted on Kay 26, 1981 (0 
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the exit temperature, T of the flowing saline-water, is low. 
With the gradual decrease in the mass flow rate and increase in 
the upward heat fluxes during the day, more amount of heat is 
to be rsnoved by less quantity of water. This results in .li^er 
exit temperature, during the later part of the day, till about 
3 p.m. After 3 p.m., further decrease in the mass flow rate is 
accompanied by a decrease in the upward heat fluxes due to low 
intensity of radiation in the evening hours. This causes a 
gradual decrease in temperature, T 

The tanperature of the saline-water inside the storage 
tank gradually rises (about 4 C - 5 C throu^out the day) due 
to absorption of heat from the surroundings, but roaches a 
steady-state after 2 p.m. 

The experimental data and the results for the second 

mode of test are given in Tables 4.2(a) and 4.2(b), respectively, 

and are plotted in Fig. 4.2. From the graph it is observed that 

the general pattern of the temperature profiles is similar to 

that in Fig. 4.1. However, the flowing saline-water outlet 

temperature, T shows a manked decrease of about 2°C - 8^C 
^0 

throughout the day. This is due to the rapid cooling of the 
water-film as its surface is being constantly exposed to the 
ambient. The temperature, ^2 of file condensation surface, is 
further reduced due to the cooling by the water-film, resulting 
in hi^er yield of distillate and better efficiency of conversion. 
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TaTole 4. 2(a) s Experimental data obtained for the Inclined 

still for the second mode of test (top ^ass 
cover removed) May 27, 1981 

© = 12 ° 



Basin 
wat er 
t emp e- 
rature 

T3_(°C) 

Glass 
c over 
t emp e~ 
ratnre 

T2(°C) 

Back 
surface 
t emp e- 
rature 

Amb i ent 
t emp e~ 
nature 

i^(°c) 

Plowing saline I 
wat er t emp eraturel 

Inlet 

I 

^i 

(^0) 

Outlet 

T 

"^0 

(°c) 

9 am 

43 

34 

34 

31 

32 

32 

10 am 

49 

40 

36 

32 

33 

34 

11 am 

54 

45 

39 

35 

34 

36 

12 noon 

58 

49 

42 

36 

35 

38 

1 pm 

60 

51 

43 ' 

37 

36 

39 

2 pm 

58 

50 

43 

37 

36 

39 

3 pm 

54 

46 

41 

37 

36 

39 

4 pm 

50 

42 

39 

36 

36 

37 

5 pm 

46 

39 

37 

35 

36 

34 
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Figures 4.3 and 4.4 show the hourly variation of still 
output obtained for the first and second modes of tests, 
respectively. From the graphs, it is found that the difference 
between the theoretical and experimental output of distillate 
for the inclined still varies from 3 pei’ cent - 20 per cent in 
case of the first mode of test and 7 per cent - 21 per cent for 
the second mode of test conducted over the entire period of the 
day, from 9 a.m. to 5 p.m. The difference in the theoretical 
and experimental results may be due to the following reasons: 


(i) loss of distillate due to leakage of water-vapour 

through the linings and dropping of condensate from 
the glass cover to the basin, 

(ii) loss of distillate during pouring and transferring 
while taking measuremaats, and 

(iii) the thin film of condensate adhering to the bottom 
surface of the glass cover has not been taken into 
account. 


Figure 4.5 shows a comparative study of the cumulative 
yield for a day for the basin type and the inclined still for 
the first mode of test. It is seen that the yields of the 


stills are nearly the same till about 12 noon, after which 


the yield from the inclined still keeps on gradually increasing 
for the rest of the day. The total increase in the yield of 
the inclined still over the basin-type unit is about 10 per cent 


From Fiff 4.6, it is found that for the second modfelcie. 

«. A VJ’H ‘ 
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when the top glass cover is removed, the total increase in the 
yield of the inclined still over the basin-type is about 38 
cent. A c caparison between Figures 4*5 aJcd 4.6 shows that the 
cumulative yield of the inclined still in case of the second 
mode of test is about 1? per cent higher as compared to that of 
the first mode. This increase in yield in the second mode of 
test is due to the facts that the bottom ^ass cover is cooled 
more effectively by the water film exposed to the ambient and 
also due to improved evaporation rate of saline water as a 
result of increased intensity of radiation falling on the basin. 

The theoretical as well as the experimental values of 
efficiency of conversion of saline water into fresh water for 
the inclined still are plotted over the day for the first mode 
of test and shown in Fig. 4.7. It is observed that the efficiency 
of conversion increases during the later part of the day. This 
may be explained as follows ; 

Since, the rate of upward heat flux (analogous to the upward 
heat losses from a collector plate to the ambient) due to the 
combined effects of convection, radiation and evaporation 
increases with the decrease in the intensity of solar radiation 
during the later part of the day, the still efficiency (q^/Ig) 
also increases. The theoretical and experimental efficiencies 
for the second mode of test are plotted in Fig, 4*8, The trend, 
in this case, is somewhat uniform throughout the day . Ccsiiparing 
the experimental results of Figures 4.7 and 4.8, it is found 
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Table 4.3 « Comparison of daily total yield of distillate 
for the two types of stills for the month of 
May, 1981 



Daily Total 
Insolation 

(KWh/m^-day) 

■ Daily Total 
y i eld of 
Inclined still 
. for 1st mode 
.of test (top 
glass cover 
retained) 

(litres/m -day) 

— ^ , 

Daily Total 
yield of 
Inclined still 
for 2nd mode 
of test (top 
glass cover 
r em cv ed ) 

2 

(litres/m -day) 

Daily Total 
yield of 
Basin-typ e 
still 

2 

(litr®/m -day) 

May 11 

4.7 

3.95 


3.45 

12 

4. 6 

- 

4. 25 

- 

13 

- 

- 

- 

3.40 

14 

4.8 

3.85 

- 


15 

4.75 

- 

4.2 

3.7 

16 

4.4 

3.80 

- 

- 

17 

4.7 

- 

4.3 

3.7 

18 

4-6 

4.05 

- 

3.45 

19 

4.75 

- 

4.3 

3.45 

20 

- 

- 

- 

- 

21 

4.90 

4.0 

- 

- 

22 

4.80 

- 

4.7 

3.4 

23 

5.10 

4.2 

- 

- 

24 

5.25 

- 

4.9 

3. 6 

25 

- 

“ 

- 

- 

26 

5.25 

4.0 

- 

3.65 

27 

5.10 

- 

4.80 

3.5 

28 

5.30 

4.0 

- 

3.6 

29 

4.90 

- 

4.60 

3.55 

30 

5.10 

3.95 


- 


5.0 

, .I- 

4.75 

3.7 




XPERIMENTAL DATA FOR THE FIRST MODE OF TEST 
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that the efficiency in the second mode of test varies hetween 
55 per cent - 70 per cent against 35 per cent - 65 per cent in 
case of the first mode. 

A comparison of the yield of distillate for the two 
modes of tests against various values of total insolation on 
an inclined still is shown in Pig. 4.9. Such a graph is useful 
in the prediction of distillate yield on an hourly, daily or 
monthly basis, for various values of the intensity of radiation. 
However, extensive data are required for accurate prediction of 
such results. 

Finally, the experimental results obtained for the daily 
total insolation and yield of distillate for the two types of 
stills during the month of May, 1981 are given in Table 4.5 and 
plotted in Pig. 4.10. 

4.3 CGHCLUSIOHS 

The following conclusions may be drawn on the basis of 
the study carried out on both the inclined and basin-type 
solar stills ; 

(1) The inclined solar still combines the process of desali- 

nation and hot-water production. The heat which is normally 
dissipated to the ambient in a conventional solar still 
here effectively to obtain hot -water on a 
at tonperatures varying between 40 °C 
throu^out the day. 
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(2) The process of cooling the bottoai glass cover by the 
flowing saline water results in high yield of distillate 
and consequently, hi^ conversion efficiency. Efficiencies 
more than 70 per cent have been attained for the second 
mode of test during the month of May. 

(3) The experimental results obtained for the inclined solar 
still match well in most cases with those obtained frcm 
theoretical analysis. The theoretical model may be 
further improved by taking into account the transient 
effects occuring inside the still. 

(4) The inclined solar still may be used for either mode of 
operation, depending upon the need of the user. The first 
mode of operation, which necessiates the use of an extra 
^ass cover is justified in view of the hot-water produc- 
tion, since it results in only a marginal increase in the 
capital cost of the still. 

(5) The major accessory of the inclined solar still is the 
saline-water storage tank, which sli^tly increases the 
capital cost. However, in actual practice, large basins 
may be constructed along sea-coasts which will be filled 
up during high tides by sea water. The saline water will 

throu^ pipes continuously into the inclined 
stills installed over avast area. Such a schane, if 
implanent ed, may result in the production of large 
quantities of freshwater economically in future. 
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AFFMIZ - 1 

EViOLUATICZ OF THE BUOYi^lJCY FACTORS FOR MEAT AIFD MASS 
TRANSFERS IE NATUB.AI, CORY EOT IGM 

Let "be the vol’ome of the enclosure occupied by the mixture 
in the inclined still. 

and , 

tn/ the total mass of the mixture 

- masses of water-vapour am sAr in the mixture, respect iv el y 
M^jMb - molecular wei^ts of wat er-va,p our and air, respectively 
p - total pressure of the mixture 

p^jPg - partial pressures of water-vapour end air, respectively 

T - mean temperature of the mixture 

densities of minture, wat er-va:o our and air, 
respectively 

Specific volumes of mixture, wat ei*— vap our , and air, 
respectively 

The coefficient of volume expansion for the mixture is given by 

e - I ( 'it ) 

Pt - -r ^ p 

Linearizing the above expression, 


P. 


1. ^ 
ir^ ^ T, - T , ^ ’ 


2 ^1 - ^2 

■srtiere, the subscri^ 1 and 2 refer to the lower and upper 
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surfaces of the enclosure, respectively, 
or, 


“^1 - "^2 


ir. 


1 2 ^ 

/i Z’ 2 

1 

-L 


/*2 

• — ^ 1 

/^1 


(i) 


Applying the ideal gas law to water-yapour and air, we obtain, 

Pa*^ = “ (ii) 


jr- V 


PbZ^= 


^ V 


where = Universal Gas Constant. 


From 

(ii), 

Pa “a 




» ir ■ 

R^T 



Prom 






ir 

V 



rrv, +•^o^n 

Adding, — = 

IT 

1 r 
E^T ^ 

Pa^a + 

Pb^b 


yf'- 1 




i. e. 

ir~ V 

t Pa“a 

+ P;^S 

] 

or , 

S7 

[ Pa^a 

+ PgMg 

] 


(iii) 
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Hence, Eqn. (i) can be written as 




‘"‘a + S’B. “b 

d d. 


“a 


Pw 

" 


] - 1 (iv) 


Erom Dalton's law of partial pressures. 


Pa + Pb 


Eqn. (iv) can be expressed as 


P.Cli-lj) =, [ 


^Ag "^A 

pT M. + (p-p. }M. 


]- 1 


A 


A 3 _" B 


[ 


Pa/a + Ag’ V^ a/a-Pa/a-Pa/b-^Pai^b 




'] 


After rearranging terms and simplify ingj 


T, 'Pa, - Pa,> <«b - “P 

P*(l,-Tg) = Tpi [1 + i = 

(ilgP - (Mg - l-ip 


] - 1 


_l r 1 + 

12 '' Mgp 


PAi - PAg 


■J - 1 


M^-M 


B A 


“ P 


An 


T -f 
12 

To 


+ [ 


^^1,^ ] fl 

^2 


3L > 


(v) 
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The coefficient of densification, is defined as follows; 


m 


_ i ( ^) 

/» '•io ’t 


(vi) 


where c denotes' the mixture concentration. The negative sign 
indicates that with the increase in the concentration of 
water-vapour in air, the mixture density, decreases, ox^ing 
to the fact that moist air is lighter than dry air. 


Linearizing the above expression, 


m 


■ ^ ) 

c. - Cp 


or. 


Pm ^"=1 " 1 




(vii) 


Appl3''ing ideal gas law to both water-vapour and air and 
proceeding in a similax may was before, we obtain, finally 


8 (Ct-c_) 

•^m 1 2 




T 


1 




(viii) 


- P 


■1 


It has been found that for small differences in temperatures 
of the two surfaces, buoyancy factors for the heat and mass 
transfers den ot ed by the left-hand sides of expressions (v) 
and (viii) give nearly the same value. 
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CO^IPUTER PRCGRAi-H-iS 


c 

c this programme computes the mass transfer rate and 

S THE HEAT FLUXES DUE TO CONVECTION , RADIATION AND 

C 

C evaporation for the inclined , SOLAR STILL C 

C C 

Qli^ttHlli:^^4^^^3^t****it***t’^*^********^****^*^ *********** **********’>^***^^^•^*'0, 

c 

c 

COMr-iPW/AV/RH0(90) ,XMU(Pa) ,XK(90) ,P(90I ,HFG(9C} ,TC90) 

CUMHOf. /BV/ A VRHO , A VMU , A VK 
REAL Nl, L , MA , »& 

DO 10 I = 27,77,1 

READ * , T(i),RHO(I),XMUfI),XK(I),P(I),riFG(I) 


C 

C 

C 

C 

C 

c 

C • 

c 

G- 

C 

C 

10 


T REPRESENTS THE TEMPERATURE IN DEGREES CENTIGRADE, HHQ, 
THE DENSITY IN KILOGRAMS PER CUBIC METRES , MU ,THS 
viscosity in kilogram PER METfiR PER SECOND , K, THE 
thermal CONDUCTIVITY IN WATTS PER METER PER DEGREE CENTI- 
CENTTGRADE , P , THE PRESSURE IN KILOGRAM FORCE PER SQUARE 
CENTTMFTEH ,AwD HFG , THE LATENT HEAT OF COMDENSATlD?< IN 
kilocalories per KILOGRAM. 

XHIUI) = XMU(I)*1.0E-05 

CONTINUE 

DO 30 N = 27,76 

DO 30 M = H+1 ,77 

M.vj s ^l ; NN s N 

CALL AVEGCMM,ijN) 

PRIi^T 20 ,T(M),TCN),AVRH0,AVMU,AVK,P(M),P(N),HFG(N) 


C 

C 

C 

C 

c 


AVRmO.AVMU,AVK REPRESENT THE DENSITY , VISCOSITY AND THERMAL 
CONDUCTIVITY OF THE AIR IN THE STILL , EVALUATED AT THE 
AVERAGE TEMPERATURE OF IRE EVAPORATING AND CONDENSING 
SURFACES , RESPECTIVELY. 

CP = 1004.10 . 


C 

C 

c 

c 


CP RFPPESEUrS THE SPECIFIC HEAT- OF THE AIR EXPRESSED IN 
JOULES PER KILOGRAM PER DEGREE CENTIGRADE, 


PRINT 25 
DO 40 TBl : 
IB « IBl-1 


1 ,41 


IB REPPESENTS THE ANGLE OF TILT OF THE SOLAR STILL NIXH Ti-i? 
HORIZONTAL , EXPRESSED IH DEGREES, 


aaoaci' 
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30 

30 

UO 

30 

00 

uO 

00 

oO 

00 

ao 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

00 

8? 

■00 

iOO 

*00 

'00 

!O0 

m 

foo 


c 

c 



t'A, ;B REPKFSEOT TriE MOLECULAR i-JElGHTS OF WATER AND ATR , 
PESPFCTIVEf.Y , PT, THE TOTAL PRESSURE IwSTDE THi:: STILL 
E^TKESSEr- IW KIuDGRA?-* FORCE SQUARE CENTIMETER. 


BbDYT = (TC 4)-T(f‘) )/TtN) + ( <P(H)-P(M) )/C(PT»MR/CMB-Mft) )-p(M J) 

1 {?(^)/ T {' 0 ) 

BULYM s (T f .n-T(F) )/T(M)4.( (p(.4)-P(H) )/(CPT»MR/(MR-MA) )-PCN)) 3=^= 
1 (T(N)/T(?1)) 


PIJOYT & BUOYM represent iE BUFYAWCY FACTORS FOR HEAT S. ^- -,5 

thahsfers , respectively, 

L = 0 , 0 b 


L REPRESENTS THE t)T STANCE IM METRES , BETWEEfJ 
AND CUN'DENSING SURFACES , 0 , THE ACCELERATION 
Hi METRES PER SQUARE SECONf 


the evaporating: 

DUE TO gravity 


G a 9 , PI 

GH a G*!.*LML*AVPH0»AVRH0M8U0YT/(AVM1JMAVMU) 

GPM a GvLML^L^AVKHOMAVRiiOfBUnYH/CAVMUMAVMiJ) 

GR.GPH REPRESENTS THE GRASHOF NUMBERS FOR HEAT AND MASS 
TRANSFERS , RESPECTIVELY 


RAT a 1 1300 ♦ (1.0 + 0.204*SIMD('4,5»(IB-37.8))) 
PR a AVMU+CP/AVK 


PR REPRESENTS THE PKANDIL NUMBER 


D a l,l7E-09M{(TCM)+TCN))/2.0+273)>lt»1.75 

D REPRESENTS THE DIFFUSION COEFFICIENT « EXPRESSED IN SQUARE 
METER PEP SECOND 


SC = CAVMU/AVRHO)/p : 

SC REPRESENT THE SCMIDT NUMBER I 

HC « AVK / LMC 2 , 44 +( 2 . 23 - 0 , 123 MIB + 0 , 00034 * IB * IB )*(( GR * PR / RAr » i 

1 COSO{FLOAT(IB)))M>Ki, 0 / 3 . 0 >)* 1 . 0 ) I 

HC REPRESENTS THE HEAT TRANSFER COEFFICIENT EXPRESSEO IN BATTS 

PER SQUARE METER PER DEGREE CENTIGRADE 


CON a HC*(T(M)-T(N3) 
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12700 
■12900 
12000 
1.3000 
13100 
1 3200 
13300 

13 400 
1 o500 
13600 
1 3700 
1 iBOO 
liOoO 
14000 
14100 
1 42 UO 
14300 

14 400 
14500 
14600 
14700 
1 49u0 
14900 
I'SOOO 
15100 
15200 
15300 

15 400 
15500 
15600 
15700 
15B00 
15700 
IbOOO 
16100 
16200 
1 6300 

1 6 4 0 0 

1 6500 
16600 
lo700 
16800 
16700 
17000 
17100 
17200 
17300 

17 4 00 
17500 
17600 
17700 
17800 
17700 
IHOOO 
16100 
1 8200 
16300 
iH400 
ins 00 
18600 
U700 

ia«oo 

18 900 
19000 



COW PfcPRFSFNTS the heat PLOX DUE TO COl'JVECTICl’^ , EXPRESSED 
WATTS PER SQUARE MEl'E*-. , 


RAO s 5.103 ♦ C((TCM)4273)/100)*»4 - ( (TfN)+273)/lOO) t* 4 ) 


RAD PEPRE5Et4TS THE HEAT FLUX DUE TO RADIATION , EXPRESSED Iw 
WATTS PER SQUARE METER 


1 CUSD(EL 0 Ar(iB)))**( 1 . 0 /J. 0 ))-l.O) 


HD REPRESENTS THE MASS -TRANSFER COEFFICIENT , EXPRESSED IN 
kilogram per square METER PER SECOND,- 


C 5 * 6 B 9 E 08 *HD/AyRHO)l=Cl./((T(M 7 4 TCN))/ 2 . 0 + 273 ))tCP(M)-pa-il 
1 (2.0=I=PT-P(M)-P(N) ) 


N1 REPRESENTS THE MASS TRANSFER RATE# EXPRESSED IN MILLILITRES 
PER HOUR 


EVAP = 0,0032283 * N1 I' HFGCN) 


EVAP REPRESENTS THE HEAT FLUX DUE TO EVAPORATION , EXPRESSED 
IN WATTS PER SQUARE METER 


PRXNT50,IB#8UOYT#BU0yM,RAr#D,GR,GRM,PR,SC,HC,HD#Ni #CnN#RA0,£-/ 

cunttnue 

CONTINUE 


F0RMAT(//132C1H=)//3X, 'fW = ' # F4, i , 4X . 'TG = ' , F4 , 1 , 4X . * RHO 

FORMATpX,I2.3X,F4.2#3X,F4.2,3X,F7.l.3X,E8.3,3X,E8,3,3X,S8,3. J 
1F4,2,3X#F4.2,3X,F4.1 ,3X,F7,5,3X#F6.1,3X»F5.i #3X,F5.t#3Xf FS. • j 
FUR.MAT(' B BUOYT BUOIM RAT • D GR GRM PH ' 


SC HC 


EVAP') 


STOP 

END 


SUBROUTINE AVEG(M#N) 

CQMMOiV/AV/RH0(90) #XMU(90) # XKC90) #PC90) ,HFG(90)vT(90) 

COMMnN/BV/AVRHO#AVMU,AVK 

AVPHO = CRHO(M) + RHOCN))/2.0 

AVMU = (XMU(M) + XMUCN)3/2.0 

AVK = CXK(M) + XK(N) )/2.0 

TXPE I 
RETURN' 

form AT OUT*) • 

END • 
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